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Abstract 

Motivated by the OPERA anomaly, we propose a matter effect enhancement of oscillation be- 
tween muon and sterile neutrinos through a neutral v s coupling to fermionic matter. We demon- 
strate that the resulting matter effect model can evade the current MINOS bound of P v „-> Vs < 0.22 
(90% CL) by participating in neutral-current interactions. We find excellent agreement between 
our sterile neutrino model parameters and parameters which fit LSND and MiniBooNE data. It is 
shown that superluminal (or " superluminal-like" ) sterile neutrinos are a viable candidate to explain 
the OPERA anomaly. 



'Electronic address: 



bramant e@hawaii . edu 



A recent result from the OPERA experiment p] measured the difference between the 
speed of light in vacuum and muon neutrino velocity along a 730 km baseline between the 
CERN CNGS beam and detection at LNGS: 

- c)/c = (2.48 ± 0.28 stat ± 0.30 sys ) x 1(T 5 [OPERA]. (1) 

This has prompted proposals of new physics[2], novel applications of standard physics to 
explain away the observed superluminal propagation [3J, and phenomenological constraints 
on violation of Lorentz invariance jH (6J [7] . 

A few proposals have utilized sterile neutrinos to model the OPERA anomaly. One pop- 
ular mechanism involves v s transport through a higher dimensional bulk [STUD] , an idea first 
promulgated in (TTJ [12] wherein active neutrinos confined to a D3 brane oscillate to sterile 
neutrinos, whose lack of gauge charge leaves them free to travel through large extra dimen- 
sions. Despite the appeal of these propositions, a striking constraint on sterile neutrino 
models (and most other superluminal neutrino models being formulated) comes from mea- 
surements of neutrinos and photons arriving from SN1987a. The detection of 24 neutrino 
events at three sites [13] arriving ~ 4 hours before SN1987a photons puts a rather stringent 
bound on superluminal electron neutrinos 

(v Ue -c)/c ~ 3xi(T 9 [IMB,KII,Baksan]. (2) 

While it might be tempting to emphasize that OPERA detected muon neutrinos and 
SN1987a produced electron neutrino data, after additional experimental constraints on neu- 
trino mass eigenstate velocity differences are applied [TUJ E] 

K-^)/c<2xl0- 19 , (3) 

active flavour-dependent velocity anomalies are ruled out (short of replacing the standard 
PMNS matrix with a different formalism [5]). 

A distinction between ([I]) and (|2j) that does survive current bounds is the energy range 
of neutrinos being detected: 

E„„ -10-40 GeV [OPERA] 
E Uc ~ 7 - 39 MeV [1MB, KII, Baksan]. (4) 

This spread in energy has motivated energy-dependent explanations of the OPERA anomaly 
in a number of papers, e.g., (TO J constructs Lagrangians with Lorentz- violating terms coming 
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FIG. 1: v s elastic scattering through a vector singlet B. 

from couplings to 10 GeV sterile neutrinos, and [H] suggests a PPW resonance [TTJ to create 
a Vf •f-)- u s oscillation sharply peaked at OPERA energies. 

However, another delineation between the SN1987a and OPERA neutrinos is that 
OPERA neutrinos propagated through a background of standard model fermions, while 
SN1987a neutrinos moved through vacuum. Motivated by these different backgrounds, in 
this paper we develop the phenomenology of a matter-dependent increase of sterile neu- 
trino production through a new v s coupling to fermions. Thus, instead of generating sterile 
neutrino mixing angles at prescribed energies through a modified energy dependence in the 
flavour mixing Hamiltonian, we instead generate a dramatic enhancement of an otherwise 
small muon-sterile mixing term in matter through a sterile neutrino coupling to standard 
model fermions. 

The contribution of the process in figure [I] to the effective potential of sterile neutrino 
propagation is given by 

</ = -S[^7^][/^T,/^], (5) 

where / is a fermion abundant in matter, namely (e~,u, d), and B is a neutral vector boson 
singlet. For the purposes of deriving a simple sterile neutrino matter effect enhancement, 
we employ this model throughout our discussion. Finally, we point out that B couples to 
(u, d)L, (v e , e)L, ur, d^, and eR with strength g$ and couples to z/ s with strength g s . The 
couplings to sterile neutrinos and fermions are assumed to be unequal owing to collider 
experiment constraints on the strength of gj (as we will see, the scale of the couplings likely 
exceeds that of electroweak interactions). 

In principle the proposed mixing enhancement will bolster all active-sterile neutrino os- 
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cillations in matter, but to simplify our analysis we assume the vacuum mixing v e>T v s is 
vanishingly small even after matter effect enhancement, so that the only appreciable sterile 
neutrino production in matter will come from •<-)■ v s . For this simple system with only 
mixing with v Sl the flavor evolution equation in matter is [15] 
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N = n e + n u + rid is the number density of matter fermions, 6 is the vacuum mixing angle 
and Am 2 is the usual squared mass difference between the mass eigenstates of z/ M and v s in 
vacuum. Note that standard model MSW terms in (|6]) have a small effect over a ~ 730 km 
baseline and have been omitted. Diagonalizing the evolution Hamiltonian yields 
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where the mixing angle and squared mass difference in matter are 
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and the corresponding matter z/„ — > v s transition probability over a distance D is 
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Before we parametrize (10) to fit superluminal propagation at OPERA, we present a 



review of applicable constraints on superluminal sterile neutrinos. In [6] it was shown that 
there is a minimum fraction of neutrinos which must travel superluminally in order to re- 
produce the OPERA anomaly. The spectral flatness of time-binned neutrino events requires 
the superluminal fraction x — ^u c+ \U^ Uc+ | 4 /£j|[/ M j| 4 to be at least ~ 0.18 at 3cr and 0.28 



at 2cr confidence. Furthermore, Cohen and Glashow [7] showed that superluminal active 
neutrinos would undergo Vf — > Vf e + e~ Cherenkov-like radiation forcing an effective energy 
cutoff above ~ 12.5 GeV for active neutrinos travelling 730 km at 7.5 km/s faster than light. 
To avoid the Cohen-Glashow energy cutoff, we might decide to require that all superluminal 
propagation occur through sterile neutrinos (as we will see below, this may not be a sufficient 
condition to avoid C-G). This requirement combined with the neutrino fraction constraint 
demands P v ^ Vs > 0.18, though a more promising model would allow for P v .„ s > 0.30. 

The most stringent bound on muon neutrino oscillation to sterile neutrinos in matter 
comes from a MINOS measurement of neutral-current (NC) interactions of the NuMI muon 
neutrino beam at the end of a 730 km baseline [16]. At face value, this MINOS result of 
802 NC events against an expected 754 ± 28 sta t ± 37 sys event background excludes P v ^ Va > 
0.22 at 90% confidence. However, in the particular case of a sterile neutrino matter effect 
there is an additional contribution to NC events from the coupling of v s to quarks ([5]). 
In fact, it is worth emphasizing that if G s ~ Gp, sterile neutrino production could both 
account for muon neutrino disappearance in prior MINOS studies [17], while mimicking 
active neutrino neutral-current interactions. As long as the mass of the intermediary boson 
B° is much greater than the momentum of the sterile neutrino, the four-fermi cross-section 
approximation <7NCu active = Gf%E u Nf should apply to sterile neutrino NC interactions as 
well. Therefore to construct matter effect models that comfortably allow P v -^ Va > 0.30, 
it is sufficient to stipulate that G s ~ Gp and m# >> E u . This stipulation has important 
implications: although we have specified that the mixing v e — > v s is very small, this is not 
well motivated, and hence G s ~ Gf reintroduces the Cohen-Glashow problem if the sterile 
neutrino oscillates to electron neutrinos. One possible conclusion is that sterile neutrinos 
are not superluminal, and some other effect accounts for the OPERA anomaly. 

3+1 sterile neutrino phenomenology has undergone development for over a decade [T8] . 
spurred by a succession of v and v short baseline anomalies and cosmological fits which 
indicate a 4th neutrino with a mass of about 0.5 eV. While this paper was in preparation, a 
conference proceedings was posted which used a similar v s matter effect to fit a 3+1 model 
to the LSND and MiniBooNE datasets [19]. Remarkably, the parameterization of [19] is 
consistent with a parameterization which allows substantial — > v s mixing at long baseline 
experiments, thereby allowing superluminal sterile neutrinos at OPERA, as we will now 
show. 
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FIG. 2: We have plotted a band of black dashes which includes the region P u ^ Us > 0.18 to 
illustrate the resonance Am 2 ~ 2EA S as detailed in the text. The band centers on the maximum 
resonance (0.975)Am 2 = 2E U A S , which is also the equation of the dashed red MINOS line (with 
E v = 3 GeV) The dotted horizontal line indicates the smallest possible Am 2 value which yields 
Py^v s ^ 0.30 at OPERA when sin 2 26 = 0.05. The dotted vertical line marks the center value 
of the sterile neutrino matter potential MINOS NC mimicry condition (G s = Gf). Underlaid is a 
scatter plot taken from [19], which fits Am^ and A s = ^/2G S N to LSND and MiniBooNE neutrino 
oscillation data. Note that the vertical axis of these dotted and dashed lines corresponds to Am 2 
of the simplified model in this paper while the scatter plot's vertical axis is Am 2 ! for a full 3+1 
model. 



In summary, a matter effect superluminal sterile neutrino explanation of OPERA requires 
Pv^v a ^ 0.30 and that the mixing angle between active and sterile neutrinos in vacuum 
is small. To enforce the latter requirement, we set sin 2 26 = 0.05, which over very long 
distances implies a vacuum z/ M — > v s transition of 0.025. Inspection of the first term in 



sin 2 IB 



(cos 20- 2V fJ G ° ) + sin 2 29 



(11) 
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produces a squared-mass difference-coupling resonance condition 

2V2NEG S 2EA S 

Am* = ~Am? (12) 

If 2 £^t >> 1, P v ^, Va will diminish rapidly, and if ^4f << 1, P v -+ Va cannot exceed a value 

of 0.05. Inserting the OPERA values into the second term of P v -> Vs , D = 730 km and E 

17 GeV, 



sin 2 
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(13) 



a minimum value of Am 2 (0.04 eV 2 ) becomes apparent. This bound follows from the maxi- 
mum value of the first P v _> Ug term (which is unity), 



2V2NEG S \ 2 sin 2 20 , . 

cos2fl- Q ~0^ 5 1 14 

Am 2 / ( nn 1-J2NKG„\ , -.2 



+ sin 

=^ P^^ s = (1) sin 2 ((0.04 eV 2 )UeV~ 2 Vom) ~ 0.30. (15) 

As Am 2 increases substantially from this value, the second term in P v ^ Vs will average to \ 
and the first term will have to resonate at ~ | to produce P„ -t Va = 0.30. 

As illustrated by the dotted horizontal line in figure (|2]), the MINOS neutral-current event 
mimicry condition (Gp ~ G s ) implies a mass-squared difference of Am 2 ~ OllO^^eV 2 . 
While at first blush it may seem odd that A s (Gp) ^ Af(Gf) ~ 10 -13 eV, the factor of 
~ 7 difference arises from the singlet vector coupling of v s to matter which avoids the (V-A) 
cancellations of the MSW matter potential. It should also be noted when examining figure [2] 
that the choice of sin 2 20 significantly impacts the location and shape of the resonance lines. 
Increasing the vacuum sterile mixing angle from the conservative value of sin 2 20 = 0.05 will 
both elevate and broaden the OPERA P v -^ Va > 0.18 inclusion band. 

If we accept the hypothesis of sterile neutrinos interacting via neutral-currents with all 
matter fermions, in the case of a straightforward singlet vector coupling the parameter space 
around the upper right portion of the MINOS resonance line is ruled out. (MINOS was not 
bombarded with NC events in [16j). However, there are intermediate values of A s along and 
beside the MINOS line, where the cross- section for sterile neutrino NC events exceeds that 
of active neutrinos, and sterile neutrinos are produced in significant quantities. A possible 
signal of this at MINOS would be an oscillation above the SM background of NC event 
counts with respect to energy (see figure [3]). 
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FIG. 3: Sterile neutrino transition probability plotted against neutrino energy for the parametriza- 
tion Am 2 = 0.45 eV 2 , A s = 1CT 11 eV, sin 2 26 = 0.05, D = 730 km. If the sterile neutrinos have a 
neutral-current interaction cross-section ~ 5 larger than SM neutrinos, an oscillation around 2.5% 
would be observable. This particular parametrization may be ruled out by virtue of its diminution 
of charged-current interactions at OPERA in the 20-30 GeV range. 

Before concluding, we observe that although the Winter j6] analysis, which fit the superlu- 
minal neutrino fraction \ — ^v c+ \U ^ Vc+ \ A /T,^!! to the flatness of the OPERA events and 
proton waveform, can be interpreted as a fraction of faster-than-light neutrinos detected by 
OPERA, it equally implies that if some new interaction is yielding muon-neutrino-like events 
and is systematically mischaracterized by St ~ 100 ns, such an interaction could account 
for the OPERA anomaly. We speculate without any immediate substantiation that owing 
to the fact that semi-sterile neutrino couplings have not been carefully studied, semi-sterile 
neutral-current interactions of sterile neutrinos could in principle mimic charged-current 
events. 

In conclusion, we have developed a matter effect enhanced model of sterile neutrino 
production at long baselines. It deserves emphatic announcement that the intersection in 
parameter space between our squared mass-coupling resonance lines and the fit of discrepant 
data from LSND and MiniBooNE, posted while this paper was in preparation, was not 
a consequence of post-hoc fine-tuning. In fact, the intersection is necessitated both by 
the neutral-current MINOS mimicry condition, without which we cannot reasonably have 
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Py^v s > 0.30, and moreover by the neutrino energies, baseline lengths, and assumption 
of relatively small vacuum P v . Vs mixing (this last assumption actually has a surprisingly 
small effect in parameter space). These necessary constraints correspond to the dotted 
horizontal and vertical lines in figure |2j There is further work on exact fits to LBL data and 
a complete 3+1 model for LBL sterile neutrino matter effect enhancement that remains to 
be undertaken. 
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